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ABSTRACT 

Context. Shock modelling predicts an electron density enhancement within the magnetic precursor of C-shocks. Previous observations 
of SiO, H"CO^, HN'^C and H'-^CN toward the young L1448-mm outflow showed an over-excitation of the ion fluid that was attributed 
to an electron density enhancement in the precursor 

Aims. We re-visit this interpretation and test if it still holds when we consider different source morphologies and kinetic temperatures 
for the observed molecules. To do this, we use updated collisional coefficients of HN'^^C and SiO with electrons in our excitation 
model. We also aim to give some insight on the spatial extent of the electron density enhancement around L1448-mm. 
Methods. We estimate the opacities of H"CO^ and HN'^C by observing the 7=3— >2 lines of rarer isotopologues. To model the 
excitation of the molecules, we use the large velocity gradient (LVG) approximation with updated collisional coefficients to i) re- 
analyse the observations toward the positions where the over-excitation of H'^^CO^ has previously been observed [i.e. toward L1448- 
mm at offsets (0,0) and (0,-10)], and ii) to investigate if the electron density enhancement is still required for the cases of extended 
and compact emission, and for kinetic temperatures of up to 400 K. We also report several lines of SiO, HN'-^C and H"CO* toward 
new positions around this outflow, to investigate the spatial extent of the over-excitation of the ions in L1448-mm. 
Results. From the isotopologue observations, we find that the emission of H'^^CO^ and HN"C from the precursor is optically thin 
if this emission is extended. Using the new collisional coefficients, an electron density enhancement is still needed to explain the 
excitation of H'^^CO^ for extended emission and for gas temperatures of < 400 K toward L1448-mm (0,-10), and possibly also toward 
L1448-mm (0,0). For compact emission the data cannot be fitted. We do not find any evidence for the over-excitation of the ion fluid 
toward the newly observed positions around L1448-mm. 

Conclusions. The observed line emission of SiO, H'^'CO^ and HN'-^C toward L1448-mm (0,0) and (0,-10) is consistent with an 
electron density enhancement in the precursor component, if this emission is spatially extended. This is also true for the case of 
high gas temperatures (<400K) toward the (0,-10) offset. The electron density enhancement seems to be restricted to the southern, 
redshifted lobe of the L1448-mm outflow. Interferometric images of the line emission of these molecules are needed to confirm the 
spatial extent of the over-excitation of the ions and thus, of the electron density enhancement in the magnetic precursor of L1448-mm. 

Key words. ISM: individual objects: L1448 - ISM: clouds - ISM: jets and outflows - ISM: molecules - physical processes: shock 



1. Introduction 

C-shock waves associated with molecular outflows are belived 
to develop a thin region of en hanced fractio n al ionisation known 
as the 'magnetic precursor' (lDrainel[l980l: iFlower et"al 
iFlower & Pineau des Foretsll2003h . Inside this region, the mag- 
netic field is gradually compressed, forcing the ions to stream 
through the neutral gas, accelerating, compressing and heating 
this fluid before the neutral one. The subsequent ion-neutral 
velocity decoupling leads to the sputtering of dust grains, in- 
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jecting large amounts of molecular material into the gas phase 
(e.g. Caselli et al. 1997; SchiUce et al. 1997; Jimenez-Se rra et al.l 
2008|) The electron density is predicted to be enhanced within 
the magnetic precursor by a factor of ~ 100 due to the fluores- 
cence UV pho tons generated after the colfisional excitation of 
H2 molecules dFlower et al.l 119961 : IFlower & Pineau des ForetsI 
2003.) . 

It has been proposed that a narrow (line-width ~ 0.5 km s ') 
and very low velocity component of SiO detected towards the 
very young LI 448 -mm outflow, is the signature o f the i nteraction 
of the magnetic precursor (Jimenez-Serra et al. 120041 hereafter 
JS04). The derived fractional abundance of SiO is of the order 
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10" in the velocity component of the precursor, an enhance- 
ment by a factor of 10 with respect to the upper limits measured 
in the quiescent ga s of dark clouds (<10~^^ in L 183 and LI 448; 
IZiurvs et al.1 Il989t iRequena-Torres et al.1 l2007h . Although the 
origin of this narrow SiO emission is still unclear, it has been 
suggested that the sputtering of dust grains at the precursor stage, 
is efficient enough to inject a considerable fraction of the grain 
mantles into the gas phase. As a consequence, bright and narrow 
SiO line emission is expected to arise from material whose cen- 
tral radial velocities a re very similar to that of the ambient cloud 
dJimenez-Serra et al.l l2b08, 2009). 

In correlation wi th the detection of narrow SiO emission, 
Jimenez-Serra et al. ( l2006l hereafter JS06) reported differences 
in the excitation of molecular ions, such as H'-^CO^, with re- 
spect to neutral molecules, such as HN'^C and H'-^CN, in the 
precursor component. In particular, the high-/ line emission of 
H'''CO+ is substantially brighter than that of HN'^'C toward 
those regions where the precursor has been detected. JS06 pro- 
posed that these differences, which cannot be accounted for by 
considering only molecular excitation by collisions with H2 and 
a single H2 density of few xlO^cm"^ for the molecular gas, 
could be produced by the selective excitation of molecular ions 
by collisions with electrons within the precursor. This study es- 
tablished that the over-excitation of the ions is consistent with 
an electron density enhancement by a factor of ~500 in the pre- 
cursor component toward LI 448 -mm (0,0) and (0,-10), which is 
simil ar to that predicted by C-shock mod elling at this dynamical 
time jFlower & Pineau des Foretsll2003h . 

However, the observed differences in excitation between 
the ion and the neutral molecules could be alter n atively ex- 
plained by opac it y effe cts. Indeed, Frerking et alj d 19791) and 
ICernicharo et al.l (Il984t) showed anomalies in the large-scale 
line intensity emission of the hyperfine components of HNC and 
HCN toward the molecular dark clouds TMC-1 and L134N. The 
equal intensities measured for these components could be pro- 
duced by the absorption of the emission arising from the dense 
cores by the less dense foreground material. This also applies to 
other high gas density tracers with high dipolar moment such as 
HCO^, since they also show similar behavi ors for the line inten - 
sity ratios between its isotopologue species (iLanger et alJll978l) . 

JS06 derived the electron density enhancement in the pre- 
cursor component of L1448-mm by assuming that the H'-'CO^, 
the HN'^C and the H'^CN emission had thin to moderate op- 
tical depths. Unfortunately, the lack of even rarer isotopologue 
observations toward this outflow prevented, first, to determine 
the actual optical depths of this emission, and second, to clearly 
establish the origin of the over-excitation of the ions in L1448- 
mm. 

In this paper we present new observations of J-3—>2 lines 
of the rare isotopologues HC''*0+, H'^C'^0+, H'^NC and 
H^^N'^C, toward the regions where the over-excitation of the 
ions has been reported. In addition, we present new observations 
of the Hnes of HN'^C and H'^CO+ toward extra posi- 

tions around the L1448-mm source, which we combine with pre- 
vious observations toward these positions, observed by JS04, to 
derive the spatial extent of the over-excitation of H'^^CO^ asso- 
ciated with the precursor These observations will allow to prove 
that the anomalous excitation of the ions is not due to a large- 
scale scattering effect, but to a real enhancement of the local 
density of electrons at the early stages of the interaction of very 
young C-shocks. 

We also re-visit the analysis of JS06, using the same exci- 
tation model as JS06 b ut including new colhsional data of HNC 
and SiO with electrons dFaure et al J2007HVarambhia et al.l2009l 



respectively), to test if the conclusions of JS06 still hold. The ex- 
citation model employed uses the large velocity gradient (LVG) 
approximation. We note that the L1448-mm outflow is likely the 
best object where the effects of coUisional excitation by electrons 
on the molecular excitation can be directly tested. Furthermore, 
in our re-analysis we fully explore the effects of source morphol- 
ogy and kinetic gas temperature on our results. 

In Sections |2] and [3] we present the observations carried out 
with the JCMT telescope and describe the results. In Section|4] 
we estimate the opacities of H'^CO^ and HN'^C toward the 
positions where the over-excitation of the ions has been re- 
ported, assuming both compact and extended source emission. 
In Section |5] we present the re-analysis of the LVG calculations 
of JS06 using the new collisional coefficients, considering dif- 
ferent source morphologies and temperatures up to 400 K. In 
this section we also analyse the data measured toward the new 
positions in L1448-mm to provide some insight into the spatial 
extent of the over-excitation of H'^CO^. Finally, discussion and 
conclusions are presented in Section|6l 

2. Observations 

We observed the 7=3— >2 transition of the rare isotopologues 
HC"*0+, Hi3Ci**0+ and H'^N'^C toward L1448-mm at off^sets 
(0,0) and (0,-10), and H'^NC toward L1448-mm (0,0). The co- 
ordinates of the central position of the L1448-mm source are 
(a = 03''25'"38^8 , 6 = 30°44'05".4 [J2000]). Finally, we mea- 
sured the 7=3^2 lines of H'^CO+ and HN'^C toward the offsets 
(0,20), (20,0) and (-20,0). 

These lines, with frequencies ranging from 248 to 267 GHz, 
were observed with the James Clerk Maxwell Telescope (JCMT) 
at Mauna Kea (Hawaii), using the position switched observing 
mode and a reference (off) position of (800", 0) with respect to 
the central position. The ACSIS spectrometer provided a spectral 
resolution of 3 1 kHz, which corresponds to a velocity resolution 
of ~ 0.04 km s"' for the A3 receiver The system temperatures 
ranged from 400 K to 700 K. 

The 7=3^2 hnes of H'^CO^ and HN'^C (at -260 GHz), 
observed with the IRAM 30m telescope at Pico Veleta (Spain) 
toward L1448-mm (0,0) and (0,-10), were taken from JS06. The 
spectra of the 7=l-*0 transitions of H'^CO+ and HN'^C, and 
of the 7=2^1 transition of SiO (at ~ 87 GHz) toward L1448- 
mm (0,20), (20,0) and (-20,0), were also observed with this tele- 
scope, and have been previously published in JS04. See JS06 
and JS04 for a full description of these observations. 

Table[T]summarises the beam efficiencies and sizes, for all of 
the observations analysed in this paper All the intensities were 
calibrated in units of antenna temperature (T^). 

3. Results 

The line profiles of the 7=3— >2 transitions of the rare isotopo- 
logues of HCO+ and HNC observed toward LI 448 -mm (0,0) and 
(0,-10), are shown in Figure[T] The 7=3^2 lines of H'^^CO+ and 
HN'^C from JS06 are also plotted. The observed parameters of 
this emission are given in Table |2] Where detected, the line in- 
tensities are well above the 3cr noise level. These data are used 
to calculate the upper limits on the opacity of the H'^CO^ and 
HN'^C emission in Section|4] 

The H'^CO^ molecular line emission toward L1448-mm 
(0,0) was fitted with a double gaussian line profile, whose radial 
velocity components are centred at V'£5r=4.7 km s"' (the ambi- 
ent gas), and at VLSR-5.2kms^^ (the magnetic precursor com- 
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Table 1. Summary of beam efficiencies and FWHM beam sizes for all the observations analysed in this paper. 



Lines 


Offset (") 


Frequency (GHz) 


Telescope 


V 


beam size (") 


Reference 


rl\^ \J {j-Z), rl U \j-Z), 

TTl5xTl3p/0 9\ H^^IVP*'^ 9^ 
n iM \^\J-Z,)^ 1 1 iy\^\J-z.) 




Z4o ~ ZD / 


J L^iVl 1 


n AO 
u.oy 


1 


This work 


TTi3pQ+/a_9\ HN^VH-?^ 




~ 260 


JCMT 


0.69 


21 




SiO(2-l), HN"C(l-0), 
H'3CO+(1-0) 


(0,0), (0,-10) 


~ 90 


IRAM 


0.82 


28 


JS06 


SiO(3-2) 


(0,0), (0,-10) 


~ 130 


TRAM 


0.74 


19 


JS06 


SiO(5-4), HN"C(3-2), 
H'^'CO+(3-2) 


(0,0), (0,-10) 


217 - 260 


IRAM 


0.52 


11 


JS06 


HN"C(4-3), H"CO+(4-3) 


(0,0), (0,-10) 


~ 350 


JCMT 


0.63 


14 


JS06 


H'^'CO+(1-0), HN'^'C(l-O) 


(0,20), (20,0), (-20,0) 


~ 87 


IRAM 


0.82 


28 


JS04 



Table 2. Observed parameters of the 7=3^2 lines of Hi3CO+, HC'**0+, H'3c"^0+, HN'^C, H'^NC and H'^N'^C toward L1448- 
mm (0,0) and (0,-10). 



(0,0) (0,-10) 





Vlsr 


Av 


n 


Vlsr 


Av 


n 


Line 


(km s-') 


(km s-') 


(K) 


(km s-') 


(kms-') 


(K) 


H'^CO+'' 


4.5(2) 


0.9(3) 


0.23(9) 


-4.7 




<0.168 




5.35(1) 


0.71(3) 


2.21(9) 


5.52(4) 


0.90(8) 


0.65(6) 


HC"*0+ 


~ 4.7 




< 0.195'' 


- 4.7 




< 0.195 




~ 5.2 




< 0.195 


- 5.2 




< 0.195 




~ 4.7 




< 0.090 


- 4.7 




< 0.087 




~ 5.2 




< 0.090 


- 5.2 




< 0.087 


HN"e 


4.75(5) 


0.6(1) 


0.35(9) 


- 4.7 




< 0.261 




~ 5.2 




< 0.237 


- 5.2 




< 0.261 




-4.7 
~ 5.2 




< 0.093 

< 0.093 








H'^N'^C 


~ 4.7 




< 0.198 


- 4.7 




< 0.207 




~ 5.2 




< 0.198 


- 5.2 




< 0.207 



" These molecular line transitions were taken from JS06. 
The upper limits of the line intensities correspond to the 3cr noise level. 



ponent). These two velocity components were already reported 
by JS04 and JS06. 

Toward L1448-mm (0,-10), the H^^CO^ 7=3->2 emission 
shows a single peak centred at the precursor component with 
VLSR=5Akms-\ The HN'^C 7=3^2 line, which is only de- 
tected toward L1448-mm (0,0), has a single peak at the ambient 
gas velocity, y£SR=4.7kms"'. The 7=3^2 lines of the rare iso- 
topologues HC"*0+, Hi3C'^0+, H'^NC and H'^N'^C were not 
detected above the 3cr limit (see Table 2). 

The profiles of the SiO 7=2— »1 lines, and of the J=3—>2 
and 1^0 transitions of H'''CO+ and HN'^C toward the L1448- 
mm (0,20), (20,0) and (-20,0) offsets, are shown in Figure |2l 
The observed parameters are given in Table[3] with all detections 
above the 3cr noise level. 

Narrow, single-peaked SiO 7=2— >1 emission is detected to- 
ward L1448-mm (-20,0), and very faintly toward L1448-mm 
(20,0). These lines are centred at the velocity of the precur- 
sor component (~ 5.2 km s '), and show line-widths of - 
0.5 km s ' (Table 3). The line profiles of H'^CO+ 7=1^0 and 
HN'^C 7=1— >0 are fitted by double Gaussian profiles, with 
line widths of 0.6 - 0.8 km s ' for each component. Toward 
L1448-mm (0,20) and (20,0), the H'^CO^ 7=1^0 lines show 
their maximum emission at the ambient velocity component with 
V£SR=4.7 km s"'. However, toward L1448 (-20,0) where the de- 
tection of narrow SiO is more evident, H'^CO^ 7=1^0 peaks at 
the precursor component with Viss=5.2 km s"' . The peak emis- 
sion of the HN'^C 7=1^0 lines are all centred at the velocity of 



the ambient gas iVLSR~4.7 km s '), except toward L1448-mm (- 
20,0) where this emission is red-shifted to -4.9 km s"' (Table 3). 

The J=3—>2 lines of H'^CO^ show single peaked line pro- 
files centred at - 5.2 km s"' toward L1448-mm (20,0) and (- 
20,0), while toward L1448-mm (0,20), a double gaussian Une 
profile can roughly be fitted. This behaviour was also reported 
by JS06 toward L1448-mm (0,0), where the 7=3^2 and 7=4^3 
lines of H'^'CO^ also show faint emission arising from the am- 
bient gas. As shown by JS06, this emission is consistent with the 
excitation of H'^CO"^ by collisions with only H2, assuming an 
H2 density of ~10^ cm"^. 

Toward (20,0) and (-20,0), the detection of H'^CO+ emission 
only in the precursor component suggests that the ions could 
have already been accelerated within the shock to larger red- 
shifted velocities. This is more apparent toward (-20,0) where 
all molecular lines have their peak emission slightly red-shifted 
to larger velocities compared to other positions. As proposed by 
JS06, this velocity shift could be understood as an evolutionary 
effect in which the molecular material toward LI 448 -mm (-20,0) 
is at a later dynamical time within the C-shock. This is also sup- 
ported by the clear detection of narrow SiO toward this position. 
We finally note that no signal above the 3cr noise level was de- 
tected for the HN'-'C 7=3— »2 lines at any position in L1448-mm 
for the precursor component. In Section 15.21 all these data will 
be analysed to establish the spatial extent of the over-excitation 
of H'-'CO^ in L1448-mm, and its implications for the magnetic 
precursor scenario. 
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Table 3. Observed parameters of the J-2-^l line of SiO, and of the 7=3->2 and /=l->0 transitions of H'^^CO+ and HN'^C toward 
L1448-mm (0,20), (20,0) and (-20,0). 







(0,20) 






(20,0) 






(-20,0) 






Vlsr 


Av 


n 


Vlsr 


Av 


n 


Vlsr 


Av 


n 


Line 


(km s"') 


(km s-') 


(K) 


(km s-') 


(km s"') 


(K) 


(km s-') 


(km s"' 


(K) 


SiO(2 ^ 1)^ 


-4.7 




< 0.027" 


~ 4.7 




< 0.030 


~ 4.7 




< 0.036 




~ 5.2 




< 0.027 


5.07(6) 


0.5(1) 


0.04(1) 


5.29(4) 


0.51(7) 


0.06(1) 


H'3CO+(l Of 


4.597(3) 


0.625(7) 


0.75(2) 


4.645(4) 


0.590(9) 


0.74(2) 


4.606(3) 


0.87(2) 


0.49(3) 




5.126(0) 


0.82(1) 


0.57(2) 


5.254(0) 


0.89(1) 


0.68(2) 


5.318(1) 


0.72(1) 


0.85(3) 


H'3CO+(3 ^ 2) 


4.5(1) 


0.5(3) 


0.20(6) 


~ 4.7 




< 0.14 


~ 4.7 




<0.09 




5.33(7) 


0.7(2) 


0.36(5) 


5.2(1) 


1.0(3) 


0.22(5) 


5.3(2) 


0.4(2) 


0.14(3) 


HN"C(1 ^ Of 


4.625(6) 


0.80(1) 


0.52(1) 


4.6(1) 


0.8(1) 


0.57(1) 


4.89(1) 


0.86(2) 


0.84(1) 




5.224(2) 


0.80(5) 


0.15(1) 


5.2(1) 


0.9(1) 


0.40(1) 


5.58(4) 


0.63(6) 


0.21(1) 


HN"C(3 2) 


~ 4.7 




<0.07 


~ 4.7 




< 0.075 


4.5(2) 


0.7(2) 


0.09(2) 




~ 5.2 




<0.07 


~ 5.2 




< 0.075 


~ 5.2 




<0.07 



These molecular line transitions were taken from JS04. 
^ The upper limits of the line intensities correspond to the 3cr noise level. 



Table 4. Isoto pic abundance ratios for the local ISM taken from 
IWilsonl([T999l) . 



Species 



Hc'^o-^ _ ric] roi 

H"cO* ["C] ["'01 

h'3c'8o+ _ 

[lOOl 



H 



H'-'CO+ 
NC _ ['^N] 



HN'-'C ["N] ["C] 
h"n"c _ ['^N] 
HN"C ["N] 



Ratio (R) 



69 
577 
1 

577 
69 
388 
1 



4. Calculation of opacities 

In order to rule out the possibihty that the over-excitation of 
H'^CO^ is due to large opacity effects, we need to prove that 
the emission lines of H'^CO^ and HN'^C reported by JS06 are 
optically thin. In addition, low optical depths associated with the 
H^^CO+ and HN'^C emission toward L1448-mm will fully jus- 
tify the use of the LVG approximation as the molecular excita- 
tion model to derive the electron density enhancement predicted 
at the precursor stage. 

We can estimate the optical depths of the H^^CO^ and 
HNi^C 7 = 3^2 lines, 732, toward L1448-mm (0, 0) and (0,- 
10) u sing the observed isotopic line intensity ratios dChoi et alJ 
119931) : 

rB.i-are 1 - e"*"^^ 



' B,32 



1 - 



This equation assumes equal excitation temperature for both 
isotopologues. Here, Tb.bi is the brightness temperature of the 
H'3CO+ 7=3^2 or HN'^C 7=3^2 Hne, and Tb,™,,. is the bright- 
ness temperature of the 7 = 3 — » 2 line of the corresponding 
rarer isotopologue (either HC"*0+ or H'''C"^0+ for H'^^CO+, or 
H'^NC or H'^N'^C for HN'^C). R is the typical isotopic abun- 
dance ratio observed in the local interstellar medium, values for 
which are given in Table |4] ( Wilson 1999). These values are av- 
erages for the local ISM. 

To estimate the brightness temperature ratios from the val- 
ues of given in Table |2] we need to know the source mor- 
phology. Since the morphology of the precursor component of 
the H'^CO+ and HN^^C emission is unknown, we assume two 
possible extreme cases: compact emission with a source size of 
5", and extended emission. 



For the source size of 5", the hne intensities in units of 
are converted to brightness temperatures by the formula Ts = 

V^^source + ^beam)/Curce whcrC Osource is the SOUrCC sizC, 0beam is 

the beam size and t] is the beam efficiency. Values for Obeam and t] 
are given in Table[T] For the extended source case, the brightness 
temperatures are approximated by T^. 

The estimated line intensity ratios and optical depths for the 
H'^CO^ and HN'^C 7=3^2 emission are given in Tables |5] 
and |6] respectively, for the different source morphologies and 
isotopologue ratios considered. 

An inspection of Tables|5]and|6]clearly shows that, many line 
ratios are larger than unity due to the lack of sensitivity of our 
observations for the rarer isotopologues. Under these conditions 
equation ([T]i cannot be used to estimate the opacities. 

From Tables |5] and |6] we draw the following conclusions: 

Ambient components: For the ambient components of the 
H'^CO^ 7=3^2 emission, the data is not sensitive enough 
to provide useful constraints. For HN'^C J=3—>2 toward 
(0,0), assuming an extended source, the ambient compo- 
nent is relatively optically thin (T32 < 1.1, derived from the 
H'^NC/HN'-'C ratio). Since the line emission from the am- 
bient component has indeed been obser ved to be extended 
dBachiller et al.lll990t ICuriel et al.|[l999l) . from our data on 
HN'-'C 7=3— >2 we conclude that the emission from this 
component is optically thin, as assumed by JS06. The am- 
bient component will not be further considered in this paper 

H'^'CO^ 7=3^2 precursor component toward (0,0): If 

the source size is > 5", this line is relatively optically 
thin (T32 < 2.1, using the estimate derived from the 
HC'**0^/H'^C0^ ratio). For the extended source case, the 
HC'**0+/H'^CO+ ratio has an upper limit of 0.09, which 
results in a negative solution for the optical depth for the 
assumed value of R. This is because equation ([T]i only has 
a positive solution for T32 for R < Ts,rare/TB,n ^ 1- Even 
if we consider the lower limit of the H'^CO^ 7=3^2 line 
intensity, corresponding to the Icr noise level, the upper 
limit of Ts.raie/TB.32 remains less than R. This suggests that 
the actual isotopic ratios in L1448 deviate from the averaged 
values of R - 69/577 in Table |4] Significant variations in 
the isotopic ratios in nearby dark clouds have been recently 
found. Indeed, the '^C/'^C ratio has been es timated to vary 
from ~ 40 to ~ 90 within the local ISM (ICasassus et al.l 
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(0,-10) 



2 
1 

(K) 

0.4 


-0.4 

0.1 



-0.1 

0.4 
0.2 


-0.2 

0.1 



-0.1 

0.2 



-0.2 




(0,20) 

I ' ' ' I 
Si0(2^1 



(20,0) 



(-20,0) 



Vlsr (km s ^) 



Vlsr (km s ) 



Fig. 1. Observations of the 7=3— >2 line of several isotopologue 
species of HCO+ and HNC toward the L1448-mm (0,0) and (0,- 
10) positions. Offsets, in arcseconds, are shown in the upper part 
of the columns. The vertical dotted lines show the ambient gas 
centred at 4.7 km s"', and the precursor component at 5.2 and 
5.4km s"' toward L1448-mm (0,0) and (0,-10), respectively (see 
JS06). 



I2005h . If we assume a extreme value of '^C/'^C^ 40, R 
decreases to ~ 0.07 providing an optical depth of T32 < 0.6. 
Therefore, our data suggest a rather low '^C/'^'C ratio and/or 
a high "^0/'**0 ratio in the L1448 molecular cloud if the 
emission is extended. In summary, the emission of H'^CO^ 
J-3—>2 for the precursor toward (0,0) is optically thin even 
if it arises from a small source. 

H'^CO^ 7=3->2 precursor component toward (0,-10): If 
the source is extended, this line is also relatively optically 
thin (T32 < 2.8, derived from the HCi'^0+/Hi3C0+ ratio). 
For a source size of 5", we can only constrain T32 to < 16. 

HN'^C /=3->2 precursor component toward (0,0) and (0,-10) 
For both geometries we cannot constrain the optical depths 
due to the lack of sensitivity of our observations. However, 
for a source of size 5", we find in Section l5.1.3l that our LVG 
model cannot reproduce the observed intensities of HN'^C, 
implying that the emission should be extended. 




2 4 6 _ 
Vlse (km s ') 



2 4 6 _ 
Vlsr (km s ') 



(km s 



Fig. 2. Observations of the SiO /=2— >1 line emission and of the 
7=3^2 and 1^0 transitions of H'^CO+ and HN'^C toward the 
L1448-mm (0,20), (20,0) and (-20,0) positions. Offsets, in arc- 
seconds, are shown in the upper part of the columns. The vertical 
dotted lines show the ambient gas at 4.7 km s"' and the precursor 
component at 5.2 km s"'. 

5. Excitation differences between ion and neutral 
molecular fluids 

5.1. Revisiting tlie ion and neutral molecular excitation 

analysis toward L1448-mm: New collisional data and the 
effect of temperature and source size 

JS06 proposed that the over-excitation of the H'^CO^ molecu- 
lar ions observed in the precursor component toward L1448-mm 
(0,0) and (0,-10) could be explained by an enhanced fractional 
abundance of electrons of ~ 5 x 10"^. We need to check if this 
conclusion still holds when using the newly rel eased collisional 
coefficients of SiO and HNC with electrons (IVarambhia et al.l 
r2009;.Faureetal 2007.) . 

In addition, JS06 assumed that the kinetic temperature of 
the gas was low [21 K, estimated from NH3 observations 
(Curiel et al. 1999)], but in this section we also explore the pos- 
sibility that the gas is at higher temperatures, which is predicted 
by C-shock models. Finally, we investigate the possibility that 
the source size is compact. 

As the excitation model, we have used the LVG code of JS06, 
but using the new collisional coefficients with electrons for SiO 
and HN'^C. It should be noted that for any isotopologues, de- 
excitation rates of the main species were employed, except for 
H'^CO^. For the latter, a proper calculation was performed and 
minor differences (less than 5%) were observed with respect to 
HCO^ rate coefficients. Excitation rates were derived from the 
detailed balance principle. 

5.1 .1 . Extended emission witli low kinetic temperatures 

The estimated H2 densities and column densities required to pro- 
duce the observed line intensities are shown in Table |7] assum- 
ing that excitation is due to collisions with H2 only. As in JS06, 
where possible we have used the 7=4— >3 and J-3—>2 lines to de- 
rive the H2 densities, as these transitions are more likely to trace 
the gas processed by the precursor compared to the lower excita- 
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Table 5. Optical depths of the H'^CO^ 7=3— »2 lines toward L1448-mm (0,0) and (0,-10) derived from different isotopologue ratios 
for the ambient (amb.) and the precursor (prec.) components, and for different source morphologies. 



(0,0) 



Ratio used 


Component 




T32(H''CO+) 


7'B,rarc/rB,32 


T32(H'-'CO+) 






Source size 5" 


Extended source 


HC"*0+/H"CO+ 


Amb. 


< 2.4 




<0.9 


< 16 




Prec. 


<0.25 


< 2.1 


<0.09 


< 0.59" 


Hi3C"*0+/H'3C0+ 


Amb. 


< 1.1 




< 0.4 


< 290 




Prec. 


< 0.11 


< 67 


<0.04 


< 24 


(0,-10) 


Ratio used 


Component 


7'B,iaic/7'B,32 


T32(H'^'CO+) 


7'B,rare/rB,32 


r32(H"CO+) 


HC"*0+/H"CO+ 




Source size 5" 


Extended source 


Amb. 


< 3.1 




< 1.1 






Prec. 


<0.85 


< 16 


< 0.3 


< 2.8 


H"C"*0+/H"CO+ 


Amb. 


< 1.4 




<0.5 


< 390 




Prec. 


<0.38 


< 280 


< 0.1 


< 83 



" This value of t(H'^^CO+) was estimated assuming R=-^, rather than R=-^ (see text for explanation). 



Table 6. Optical depths of the HN'^'C 7=3— >2 line toward L1448-mm (0,0) and (0,-10) derived from different isotopologue ratios 
for the ambient (amb.) and the precursor (prec.) components, and for different source morphologies. 



(0,0) 


Ratio used 


Component 


TB.r.JTB.31 r32(HN"C) 


7'B,rarc/7'B.32 


T32(HN"C) 






Source size 5" 


Extended source 


H'^NC/HN^C 


Amb. 


< 0.74 < 7.6 


<0.3 


< 1.1 




Prec. 


~ 1.1 


-0.4 


-2.5 


H"N"C/HN"C 


Amb. 


< 1.6 


<0.6 


< 320 




Prec. 


-2.3 


-0.8 


-700 


(0,-10) 


Ratio used 


Component 


TB.r.JTB.31 r32(HN"C) 


7'B,rarc/7'B.32 


r32(HN"C) 






Source size 5" 


Extended source 


H'5N'3c/HN'3c 


Amb. 


-2.2 


-0.8 


- 610 




Prec. 


-2.2 


-0.8 


- 610 



tion lines. Where this information is not available, we have used 
the J-3—>2 and 7=1— >0 lines to estimate the H2 densities. The 
H2 densities for all molecules are a factor of ~ 10 times higher 
than those in Table 3 of JS06. This difference mainly arises be- 
cause JS06 used main beam temperatures to estimate the H2 den- 
sities, which is equal to Tb when the source size exactly fills the 
beam, whereas in our case we have used T^, appropriate for ex- 
tended emission. 

We find that an H2 density of ~ 10^ cm"^ can explain the 
excitation of the neutral species, but the H2 density required for 
H'^CO^ is a factor of ~ 5 times larger toward (0,0), and ~ 10 
times larger toward (0,-10), as previously found by JS06. 

If we now assume an H2 density of 1 x 10^ cm"^, fractional 
ionisations of 5 x lO"'* and 7 x 10"^ toward (0,0) and (0,-10), 
repectively, are required to explain the line intensity ratios mea- 
sured for H'^CO^ in the precursor component. These enhanced 
electron abundances do not make a strong impact on the pre- 
dicted intensities of SiO and HN'^C, since they differ by less 
than 30 per cent for SiO, and by a factor of 2 for HN'^C, with 
respect to the observed ones. This shows that when using the new 
coUisional coefficients of SiO and HN'^^C in the calculations, the 
electron density enhancement in the precursor component is con- 
sistent with the measured intensities of both the ion and the neu- 
tral species. The LVG model estimates the column densities and 
opacities of the H'^CO+ 7=3^2 Hnes to be = 1 x 10^^ cm"^ 
and 732 = 0.2 toward (0,0) and = 3 x 10" cm^^ and T32 = 0.05 



toward (0,-10). These opacities are consistent with the observed 
upper limits in SectionH] 

5.1.2. Extended emission witli liigli kinetic temperatures 

Since the coUisional coefficients forH'^CO^ with electrons used 
so far in the LVG model are only available for transitions up to 
7 = 5, which lies at ~ 62.5 K, for the LVG calculations at higher 
temperatures w e use updated coUisional coefficients for HCO^ 
with electrons (iFaure & Tennvsonll2001l) . with data for transi- 
tions up to 7 = 20, lying at ~ 900 K. To ensure that the LVG 
calculations are converged with respect to the rotational popu- 
lations, we can only perform LVG calculations for temperatures 
< 400 K. 

If the gas has a kinetic temperature of 100 K, the estimated 
H2 densities for the neutral species fall to ~ 3 x 10^ cm"^^. For 
H'^CO^, the estimated H2 densities are ~ 4 x 10^ cm"^ and 
~ 8 X 10^ cm"^ toward (0,0) and (0,-10) respectively. Assuming 
that the gas density is 3 x 10^ cm"^, the H'-'CO^ ions would still 
require fractional ionisations of 3 x lO"** and 8 x 10"^ toward 
(0,0) and (0,-10) respectively to fit the observations. 

Performing a similar analysis at 200 K, 300 K and 400 K, we 
find that toward (0,0), the required H2 density to fit the H'^CO+ 
data is a factor of 1 .5 -2 times larger than the H2 density required 
for the neutral species, and toward (0,-10) this factor is 2.5 - 
4. The fractional ionisations required to fit the data assuming 
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Table 7. H2 densities and molecular column densities of H'^CO^ and HN'^C derived assuming extended emission and a temperature 
of 21 K for all molecules, for the precursor component of L1448-mm. 





H2 density (cm ^) 


Column density (cm *^) 


Molecule 






(0,0) 


H"CO+ 


5.0 X 10" 


9.3 X 10" 


SiO 


1.6 X 10*^ 


1.2 X 10" 


HN"C 


< 1.3 X 10** 


3.7 X 10" 


H"CN 


< 4.8 X 10'^ 


3.3 X 10" 


(0,-10) 


H"CO+ 


3.7 X 10' 


2.6 X 10" 


SiO 


3.5 X lO** 


1.4 X 10" 


HN"C 


< 7.1 X 10^ 


6.8 X 10" 


H"CN 


< 3.6 X 10** 


4.2 X 10" 



only one H2 density, range from 3xl0^to2xl0'' toward 
(0,0), and from 8 x 10""^ to 4 x lQ-\ toward (0,-10). Without 
the electron density enhancements, the H'^CO^ 7=3— >2 lines 
are underpredicted by 20 - 60%, and the H'^CO+ 7=4^3 lines 
are underpredicted by 30 - 70%. The estimated column densities 
and opacities of H'^CO^ 7=3^2 are slightly lower than for the 
21 K case above, with = 8 - 9 x 10'' cm"^ and T32 = 0.2 
toward (0,0) and = 2 x lO" cm"^ and T32 = 0.03 toward (0,- 
10). Again these opacities are consistent with the upper limits 
derived from observations (see Section 4). 

Note that toward (0,0), even though the neutral species re- 
quire an H2 density of only 1.5-2 times less than that for 
H'^CO^, the estimated fractional ionisation is very high, up to 
four orders of magnitude higher than the typical fractional ioni- 
sation in dark clouds (~ 10~^). This is because the electrons do 
not significantly affect the excitation of H'^CO^ until they reach 
a fractional abundance of a few xlO^^*. This is demonstrated by 
Figure[3] where we show the H'^CO^ 7=4^3/7=3— >2 ratio as 
a function of the electron fractional abundance predicted by the 
LVG model, for an H2 density of 3 x lO*" cm""*, and a tempera- 
ture of 100 K. For each position we have used the column den- 
sities estimated above, although Figure [3] shows that variations 
in the column density have a minimal impact on the predicted 
7=4^3/7=3^2 ratios. 

We note that considering the Icr errors in the observa- 
tions, toward (0,0) the lower limit of the observed H'^CO^ 
7=4— »3/7=3^2 ratio is approximately equal to that predicted 
by the LVG model for low electron densities (X{e)< 10"^) for 
temperatures of 100-400 K. It is therefore possible that there is 
no need for an electron density enhancement toward (0,0) for 
these temperatures. However, toward (0,-10), as demonstrated 
by the error bars plotted in Figure |3] the need for an electron 
density enhancement to explain the observations is more com- 
pelling; the lower limit of the observed 7=4— »3/7=3^2 ratio is 
13-24 per cent larger than that predicted by the LVG model for 
low electron densities for these high temperatures [compared to 
within 0.2-12 per cent toward (0,0)]. 

5.1.3. Compact emission 

Assuming a source size of 5", converting the T^s to brightness 
temperatures, as described in Section|4] gives brightness temper- 
atures of40-50KfortheHi^^CO+ 7=1^0 lines, and 15-30K 
for the HN'^C 7=1^0 Hues. Therefore, if the emission is com- 
pact, the kinetic temperatures of the gas must be > 50 K. The 
brightness temperatures for the H'^CO^ 7=3^2 and 7=4^3 



1.5 



1.0 



0.5 



1 1 1 1 1 1 
observed ratio (0,-10) 


1 ' ^^^^^ ' 


observed ratio (0,0) t 


7 


1 ... 1 ... 1 . 



-8 -6 -4 -2 

log[X(e)] 

Fig. 3, The U^^CO^ 7=4^3/7=3-^2 ratio as predicted by the 
LVG model at 100 K toward L1448-mm (0,0) (solid curve) 
and L1448-mm (0,-10) (dotted curve), against the electron frac- 
tional abundance. The observed 7=4— >3/7=3— »2 ratios of the 
pre-cursor components towards both positions are plotted, and 
the eiTorbars coiTespond to the Icr eiTors in the observations. 



Unes are only 20 - 25 K toward (0,0) and ~ 7 K toward (0,-10). 
For the HN'^C 7=3— >2 lines, their brightness temperatures are 
< 3 K toward both positions. For gas temperatures > 50 K, our 
LVG model cannot fit the HN'^'C intensities toward (0,0) and (0,- 
10), assuming that the molecular excitation is produced by col- 
lisions with only H2. This is because we cannot fit the extremely 
low values of the 7=3— >2/7=1^0 brightness temperature ratios 
for HN'^C. In fact, the LVG model cannot reproduce the HN'^C 
intensities even for kinetic temperatures of 15 - 30 K. There is 
a similar problem for H'^CO^; whilst it is possible to find an 
H2 density and H'^CO^ column density to fit the 7=4— >3 and 
7=3— »2 observations for kinetic temperatures > 50 K, these fits 
cannot match the extremely high intensities of the 7= 1 — >0 lines, 
and it is impossible to fit the extremely low 7=3^2/7=1— >0 
line intensity ratio. Including an electron density enhancement 
increases the 7=3— >2/7=1^0, making the difference between 
the model and observed values even larger Furthermore, obser- 
vations of NH3 toward L144 8-mm indicate tha t the neutral gas 
has a temperature of ~ 20 K dCuriel et al.ll 19991) . lower than the 
minimum temperature of ~ 50 K required for compact emission. 
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The observational evidence and our modelling results are there- 
fore not consistent with compact emission. 

5.2. Analysis of new positions around the L1448-mm source 

In this section, we present the analysis of the molecular data to- 
ward the newly observed positions L1448-mm (0,20), (20,0) and 
(-20,0). We use again the LVG approximation as the molecular 
excitation model. The ion and neutral low- and high-7 lines of 
H'^CO^ and HN'^C measured toward these offsets, are needed 
to constrain not only the spatial extent of the emission arising 
from the magnetic precursor, but also the extent of the electron 
density enhancement associated with this early shock stage. As 
suggested by JS06 for the L1448-mm (0,0) and (0,- 10) positions, 
if an electron density enhancement occurs in the precursor, the 
high-7 rotational lines of H'^CO^ are expected to be signifi- 
cantly brighter than those of HN'^C. 

If the emission is optically thin (see Section |4]l, the line in- 
tensity ratio between different transitions is directly related to 
the excitation temperature of the molecular species. In Table [8] 
we show the 7=3^2/1 ^0 integrated intensity ratios for H'^CO^ 
and HN'^C toward all positions measured in L1448-mm. 

For H'^CO+, we find that the 7=3^2/1^0 ratio is signifi- 
cantly larger in the precursor component than in the ambient gas 
(by more than a factor of 3), only toward L1448-mm (0,0) and 
(0,-10). For HN'-'C, this ratio is similar (i.e. differs by less than 
a factor of 3) in both velocity components toward all positions, 
suggesting that the over-excitation of H'^CO^ does not show a 
large-scale morphology around the L1448-mm core. This clearly 
contrasts with the idea that the anomalous excitation of H'^CO^ 
is produced by a diffusion/scattering effect (Frerking et al. 1979; 
Cernicharo et al. 1984). 

This can also be qualitatively shown by deriving the H2 
densities required to explain the line intensity ratios shown in 
Table [8] as in Section 15.1.11 We assume that the emission is 
spatially extended and consider a kinetic temperature of 21 K. 
The estimated H2 densities and molecular column densities are 
shown in Tabled 

We see that the intensities of the low- and high-/ H'^CO^ 
lines arising from the precursor component toward L1448-mm 
(0,20), (20,0) and (-20,0), can be explained by an H2 density that 
differs only by a factor of 2 with respect to that of the ambient 
gas, and that is less than the upper limits of the H2 denstities 
derived from HN'^C within the precursor component. 

These results suggest that the ion H'^CO^ fluid does not 
show any evidence of over-excitation with respect to the neutral 
fluid toward the newly observed offsets around L1448-mm. This 
implies that molecular excitation by collisions only with H2, can 
explain the intensities of the low- and high-7 lines of H'^CO^ 
and HN'^C measured toward these new positions. The detected 
over-excitation of the ions toward L1448-mm (0,0) and (0,-10) 
seems to be restricted to a particular region located in the south- 
ern red-shifted lobe of this outflow. This is the ex pected behavior 
since v ery young shocks have been reported bv lGirart & AcordI 
(1200 ih toward these positions. 

6. Discussion and conclusions 

Anomalies in the intensities of the hyperfine components of 
the strongly polar HCN molecule were reported on larg e-scales 
toward the TMC-1 molecular dark cloud by Cernich aro et alJ 
(fT98l . These authors showed that these anomalies are due to 
large opacity effects and are produced by the re-emission of ra- 



diation coming from the densest part of the core, by a more ex- 
tended and less dense envelope. If this applies to the L1448- 
mm case, then we should expect to detect such differences in 
the excitation of other strongly polar molecular species such as 
HCO^ and HNC, on large-scales around the L1448-mm core. 
In Section H] from isotopologue observations we have shown 
that the precursor component of H'^CO^ is optically thin toward 
(0,0), and if the emission is extended, this component is also op- 
tically thin toward (0,-10). This implies that the over-excitation 
of H^^CO^ in the precursor toward these positions cannot be due 
to a large-scale optical depth/scattering effect. 

In view of the recent relea se of new collisional coeffici ents of 
SiO and HNC with electrons (iFaure et al. 2007; Varambhia et al.l 
2009), we have re-analysed the multi molecular line observa- 
tions toward L1448-mm (0,0) and (0,-10) of JS06, but taking 
into account different possibilities for the source size extent and 
temperature. Assuming extended emission and kinetic tempera- 
tures of up to 400 K, and considering that the molecular exci- 
tation is due to collisions with only H2, a higher H2 density is 
required to match the observations of H'^CO^ than that needed 
for the neutral species. For temperatures of 100 - 400 K, toward 
L1448-mm (0,0) it is possible that, considering the Icr errors in 
the observations, the observations can be explained by a single 
H2 density for the ions and the neutrals, but toward (0,-10) these 
errors are not enough to explain the difference in the required 
H2 densities. This leads to the need for an extra excitation mech- 
anism for H'^CO^ in the precursor component of L1448-mm. 
This extra excitation can be provided by an electron density en- 
hancement of up to a factor of ~ 10"*. As in JS06, the derived 
electron density enhancement exceeds those predicted by MHD 
shock mo delling by up to two orders of m agnitude ( Flower et al.l 
1996; Fl ower & Pineau des Foretsll2003h . However, simply con- 
sidering the Icr errors in the observations can reduce the electron 
density estimates by an order of magnitude. Due to these very 
large uncertainties, we can conclude that an electron density en- 
hancement can explain the observations, but we cannot quantify 
it with an accuracy better than one order of magnitude. 

If the emission is compact, the gas kinetic temparture must 
be > 50 K. However, we find that the HN'^C and H'^CO+ lines 
cannot be fitted in this case, so we conclude that the emitting 
region associated with the precursor must be > 5". 

The observations toward the L1448-mm (0,20), (20,0) and 
(-20,0) positions, do not show evidence for an over-excitation 
of H'^CO^ with respect to the neutral HN'^C molecules in the 
precursor component of LI 448 -mm. This suggests that the extra- 
mechanism responsible for the over-excitation of the ions has its 
origin in local phenome na likely linked to the r ecent interaction 
of very young C-shocks. Girart & Acord (2001) have indeed re- 
ported very young shocks (with a dynamical age of < 90 yr) to- 
ward the southern position L1448-mm (0,-10), making the prob- 
ability to detect the magnetic precursor of C-shocks not negli- 
gible. Since our results show that the over-excitation of ions is 
confined to this region, there seems to be strong evidence to sug- 
gest that this over-excitation is produced by the electron density 
enhancement at the magnetic precursor stage of young shocks. 
Furthermore, the confinement of the over-excitation of the ions 
also implies that this effect is unlikely to be caused by external 
UV illumination, as it should be observed tracing the inner re- 
gions of the outflow cavity, and on larger spatial scales. 

In order to fully understand the extent of the magnetic pre- 
cursor of L1448-mm, high angular resolution observations of 
the high-7 transitions of H'^CO+, HN'^C and SiO caiTied out 
with interferometers such as the Submillimeter Array (SMA), 
are therefore strongly required to measure the source sizes of this 
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Table 8. 7=3— > 2/1^0 integrated intensity ratios. Data were taken from JS06 and this paper. 





(0,20) 


(20,0) 


(-20,0) 


Molecule 


Ambient Precursor 


Ambient Precursor 


Ambient Precursor 


H"CO-" 


0.2 0.5 


< 0.2 0.4 


< 0.2 0.09 


HN"C 


<0.1 <0.3 


<0.1 <0.1 


< 0.1 < 0.3 




(0,0) 


(0,-10) 


(0,-20) 


Molecule 


Ambient Precursor 


Ambient Precursor 


Ambient Precursor 


H"CO+ 


0.5 1.6 


< 0.1 0.6 


< 0.5 0.2 


HN"C 


0.4 < 0.2 


< 0.2 < 0.1 


< 0.2 < 0.3 



Table 9. H2 volume densities and molecular column densities of H'^CO"^ and HN^^C derived assuming extended emission for both 
the ambient and the precursor components of L1448-mm. 



H2 density (cm ^) Column density (cm 



Molecule 


Ambient 


Precursor 


Ambient 


Precursor 


(0,20) 




1.4 X 10^ 


3.5 X l(f 


4.8 X W 


4.0x10" 


HN"C 


< 1.9 X 10' 


< 9.2 X 10' 


< 6.4 X 10" 


< 1.5 X 10" 


(20,0) 


H"CO-H 


< 8.8 X lO"* 


1.8 X 10' 


< 5.1 X 10" 


4.3 X 10" 


HN"C 


< 1.8 X 10' 


< 3.3 X 10' 


< 7.2 X 10" 


< 4.1 X 10" 


(-20,0) 




< 8.7 X 10* 


7.1 X 10" 


< 3.5 X 10" 


6.4 X 10" 


HN"C 


1.2 X 10' 


< 6.1 X 10' 


1.4 X 10'^ 


< 2.0 X 10" 



molecular emitting regions and to confirm whether or not this 
emission is optically thin. Observations of very high-7 (J > 10) 
transitions are also required to verify whether or not the gas has 
reached the high temperatures predicted by C-shock models, as 
these fines should be detected for such temperatures. It is also 
necessary to search for signs of the electron density enhance- 
ment toward other young outflows where narrow SiO has been 
detected, such as L1448-1RS3, NGC1333-IRS4 and NGC1333- 
IRS2, in order to investigate whether the electron density en- 
hancement is a common phenomenon in this type of objects, and 
if it is correlated with the detection of very narrow SiO emission 
centred at ambient cloud velocities. 
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